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Abstract 

Lieb, Schultz and Mattis (LSM) [T] studied the S=l/2 XXZ spin 
chain. The theorems of LSM’s paper can be applied to broader models. 

In the original LSM theorem the nonfrustrating system was assumed. 
However, reconsidering the LSM theorem, we can extend the LSM the¬ 
orem for frustrating systems. Next, several researchers tried to extend 
the LSM theorem for excited states. In the cases = ±1, ±2 • • •, the 
lowest energy eigenvalues are continuous for wave number q. But we 
found that their proofs were insufficient, and improve upon them. In 
addition, we can prove the LSM theory without the assumption of 
the discrete symmetry, which means that the LSM-type theorems are 
applicable for Dzyaloshinskii-Moriya type interactions or other non- 
symmetric models. 

Keywords: Lieb-Schultz-Mattis, rigorous theorem, frustration, one- 
dimension, Dzyalosinskii-Moriya 


1 Introduction 

In statistical physics, exact solutions such as Onsager’s theory for a two- 
dimensional (2D) Ising model, the Bethe ansatz for ID quantum systems 
p El n E], the matrix product method PE], etc, have played important 
roles. 

Besides exact solutions, there are rigorous theorems such as the Mermin- 
Wagner theorem, the Marshall-Lieb-Mattis theorem the Lieb-Schultz- 

Mattis theorem [T] etc, which do not give quantitative, but qualitative re¬ 
sults. Since such rigorous theorems are based on symmetries, they can be 
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applied to broader models. And one can use them to check the consistency 
of approximations, experiments, or numerical data. 

Lieb, Schultz and Mattis (LSM) [1] studied the S=l/2 XXZ spin chain. 
In appendix B of their paper, they described two theorems. In the hrst 
theorem, it was proved that the ground state is unique in hnite L systems. 
In the second theorem, they proved that there exists a low-energy 0{1/L) 
excited state; in the inhnite limit, this means that either there are degenerate 
ground states or a vanishing gap. 

The hrst theorem was nothing more than an extension of Marshall’s the¬ 
orem [8], and was later applied to more general cases (higher dimensions, 
bipartite lattice, ferrimagnetism, etc.) by |H], therefore it is appropriate to 
call the hrst theorem the “Marshal-Lieb-Mattis (MLM) theorem”. The MLM 
theorem was applied to the spin model with biquadratic exchange 110] . Since 
the Perron-Frobenius theorem is used in the MLM theorem, it is limited in 
the non-frustrated case. 

The second LSM theorem was extended for general spin S and was applied 
for various models by nn, and it was proven that there exists a low-energy 
0{1/L) excited state for half-odd-integer spin cases. It is worth noting that 
the assumptions of the MLM theorem and those of the second LSM theorem 
are independent, though this fact has been overlooked. Hereafter we call the 
second LSM theorem alone the LSM theorem. In this paper we shall extend 
the LSM theorem without the assumption of the uniqueness of the ground 
state, using a squeeze theorem type discussion. Therefore, we can extend the 
LSM theorem for frustrated systems. 

On the other hand, Kolb [12] studied the energy spectra of the XXZ spin 
chain with the twist boundary condition, independent of mini. He obtained 
similar conclusions regarding the two-fold degenerate ground state for the S 
half-odd-integer case. He also pointed out the change of the wave vector 
when varying the twist boundary condition for nonzero total spin cases. 
In section H of [13], combining Kolb’s idea and LSM theorem, Fath and 
Solyom insisted that they proved the continuity of the dispersion curve for 
S'! = ±1, ±2 • ■ ■. However, we hnd that their proof is insufficient. 

Oshikawa et al. [TT], using the LSM type discussion, pointed out that 
there may be magnetic plateaux. In addition, they emphasized the impor¬ 
tance of either the space inversion or the spin reversal symmetry, besides 
the U(l) and the translational symmetry. Although the proof of the LSM 
theorem becomes simplihed with the assumption of discrete symmetry, it 
excludes nonsymmetric spin models. 


2 


In this paper we prove the LSM theorem, by using only the assump¬ 
tions of the U(l) symmetry, the translational invariance and the short-range 
interaction. There is no need for the uniqueness of the ground state for a 
hnite system, the space inversion or the spin reversal symmetry assumptions, 
which means we can extend the LSM theorem for frustrating or nonsymmet- 
ric models (e.g. Dzyaloshinskii-Moriya interaction HSUS])- 

The layout of the paper as follows. In section 2, we introduce the deh- 
nition of symmetry operations. Section 3 is the main part of this work: we 
prove the continuity and the periodicity of the lowest energy spectra as a 
function of wave number q, assuming the U(l) and the translational symme¬ 
tries plus the short-range interaction. Furthermore, in section 4 we discuss 
the discrete symmetries, i.e. the space inversion and the spin reversal. We 
consider several specihc spin models with various lattice structures and sym¬ 
metries in section 5. In section 6, using the nonfrustration condition (model 
specihc nature), we can discuss the minimum location of spectra. In section 
7, we compare our theorems with DEI and [H]. In section 8, we illustrate 
possible spectra for frustrating systems, those for magnetization plateau, and 
those for nonsymmetric case. Section 9 is the conclusion. 


2 Model, symmetries, eigenstates 


In this section we consider the symmetries of the spin chain. As a typical 
model, we treat the following ID generalized XXZ spin Hamiltonian: 


L i/2 
j=l r=l 


J(r) 




+ A(r),VS: 


j j+r 


+ '‘XA 

i=i 


( 1 ) 


where {SjY = S'(S' -|- 1) (S' = 1/2,1, • • •), with the system size L even and 
the periodic boundary condition (PBC): 


QX,y,z _ QX^y^z 

^L+j - 


( 2 ) 


Hamiltonian ([T]) is invariant under the space inversion. 

And when h = 0 in ([1]), the Hamiltonian is invariant under the spin 
reversal. 
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2.1 Symmetries 

Next we enumerate the symmetry operations. Hereafter we denote 

ILL 

QX _ \ ^ QX QV _ \ ^ oy QZ _ \ ^ QZ 

- 2^ >=’j, — 2-^ 

j=l j=l j=l 

1. Rotation around the ; 2 -axis 

We dehne the rotational operator of the ; 2 -axis as 

= exp(-m^), 

C/«‘ = S* exp(±W), (U^YS’U^ = S’. 

2. Translation operator by one-site: Htri- 


( 3 ) 


(4) 

(5) 


= S'fl 


x,y,z 




( 6 ) 


3. Space inversion (site parity). 


pt _ gx,y,z 


(7) 


Because = 1, it is shown that P ^ = pl = P and the eigenvalue of 
P is ±1. 


4. Space inversion (link parity). 


pt qx,y,z p _ qx,y,z 

-'link '-’i—j-j-l- 


( 8 ) 


The link parity operator can be dehned as the product of the site parity 
operator and the translational operator. 


Plink — PPtrl- 


(9) 


5. Relation between site parity and translation. 

There is a relation between parity and the translation operation: 

Pu,P,P = U2l ( 10 ) 

therefore 

(PPtrl)' = 1, (11) 

that is, the eigenvalue of link parity is also ±1. 
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6. Spin reversal 

The operator of tt rotation around the y-axis is given as 

= exp(-7ri5f,) (12) 

then 

(UiVSfUi = -SJ, (U«JS]Ui = -S]. (13) 

The eigenvalue of the operator is an integer from the evenness of 
L. Therefore, we obtain that = 1 and 
that the eigenvalue of is ±1. 

7. The operators are invariant under the translation and the space 

inversion. In addition, is invariant under the rotation around the 
axis. 

2.2 Vacuum 

We take the fully aligned spin state as a vacuum: 

4^10) = ^|0), .S+|0) = 0, (0|0) = 1. (14) 

2.3 Eigenstates 

We write the eigenstate for the translation and the total spin as 

S^\S^-,q) = S^\S^-,q), U,,,\S^-,q)=exp{tq)\S^-,q). (15) 

Moreover, when the Hamiltonian is translational and U{1) invariant, one can 
write 

H\S^-,q) = E{S^-,q)\S^-,q) (16) 

2.4 Theorem on the translation, space inversion and 
rotation around the y axis in PBC 

We review several theorems of the PBC case. 

Theorem 1. 1. Energy spectra are 27r periodic with wave number q. 

E(^-g + 27r) = E(^-g). (17) 
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2. For the wave number g 7^ 0, vr, by using (ITU]) , one can show 


P\S^-q) = \S^;-q). (18) 

Parity of the eigenstate g = 0 or g = tt is well dehned: 

P\S^-,q = 0) = ±\S^-,q = 0), 

P\S^-q = 7r) = ±\S^-q = 7r). (19) 

3. For S'!. 7 ^ 0 , by using f[T3|) . one can show 

UyjS^-,q) = \-S^-,q). ( 20 ) 

The eigenstate of S'!. = 0 is also the eigenstate of with eigenvalue 
il. 

UyjS^ = 0-,q) = ±\S^ = 0;q). ( 21 ) 

4. 

(ppyp = ppyppy = pypy = i ( 22 ) 

thus the eigenvalue of Ppy is ± 1 . 

5. When the Hamiltonian is invariant under the space inversion, 
energy spectra are symmetric about g = 0 : 

E{Sp-q) = E{SPq). (23) 

6 . When the Hamiltonian is invariant under the spin reversal, 
energy spectra are symmetric under the spin reversal: 

E{-SPq)=E{SPq). (24) 


3 Extension of the LSM theorem 

In this section, we will extend the LSM theorem without the assumption of 
the uniqueness of the lowest energy, by using squeeze theorem type methods. 
And we will use only the translational and the U(l) symmetry, and will not 
treat the Hamiltonian directly. We do not assume the discrete symmetries 
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such as the space inversion or the spin reversal. Hereafter we express \S^; q) 
as one of the lowest energy eigenstates in the subspace of and q. 

We dehne the twisting unitary operator as 

USi, = - Syj , (25) 

then we obtain 

(C/S,)'s;C/5, = S+exp{±2mj/L), 

(USiJS-USi, = S-exp(T2ir!i/L), 

((7SjtS|C/5, = S', (26) 

and 

[US)-' = (I7‘:)t = (7“,, (27) 

and 

t/£|0) = |0). (28) 

Doing unitary transform ([1]) with twisting operator, we obtain 

- H 

ll/2 

= 5^5^ -p(>^/>^7+r(exp(T27rri/L) - 1) + h.c.). (29) 

j=l r=l 

3.1 Main theorem 

Lemma 1. (Translation operator and twisting operator) 


^ /s ^ ^ / 27r? 

= t^nit/r2.exp f ±^(^T 

^ A ^ / ‘2717 

= 17S.f^triexp f T^(^t 



(30) 
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Proof. 


PfiSiAri = exp E - ■») 

= exp (eO- - 1)(^I 

= t>5,exp('±Hp(S|.-SL) 


- S) + L(S£+, 
exp(=F27ri(^^ 



(31) 


where we used = Sf. Combining this equation with the fact that the 
eigenvalue of — S is an integer, we obtain fl5U]l . 

□ 


Theorem 2. In the subspace with a quantum number Sf, on the lowest 
energies of the three wave numbers q, q ± 2'KSf/L + 271S, the following 
inequality holds: 


E{Sf; q - 27iSf/L + 27iS) + E{Sf; q + 27iSf/L + 27iS) - 2E{Sf-, q) < 0(1/L). 

(32) 


Proof. The following combination 

+ {UXVHU%^ - 2H (33) 


is translational invariant from lemma 1. 

And from lemma 1, we obtain 

UUUT 2 .\Sf-, q)) = exp(z(g ^ + 27iS)){Uf,^\Sf- q)) (34) 

since 2S is an integer. 

Using fl33l) and fIM)) . we can prove the following inequality: 

E{Sf- q - 2nSf/L + 27r^) + E{Sf] q + 27iSf/L + 2nS) - 2E{Sf] q) 

< (Sf; + {U%JHU%, - 2i/)|SJi q) 

L 

= E E ^(>')(cos(2xr/L) - 1){S|,; 91(5+5/+, + 5/5/+,)|5|.; 9 ) 

j=l r=l 

< 0 ( 1 /L), 


(35) 





where we used the variational principle, and the next relation 

(S|.;«|(S+S-+, + S-S+,)|S|;«) = (S|.;«|(S+Sr+, + SrhV)l'!^f;^ (36) 

using translational operations, and we assume that the transverse interaction 
is short-range (for example |J(r)| oc exp(—m|r|). More detailed discussion 
on the interaction range is found in [19]. □ 

[Remark] 

• The longitudinal interaction A(r) and the magnetic held h give no 
restriction on theorem 2. 

• Although the form in line 3 of (l3^ seems specihc for the model ([!]), 
one can show (l3^ for multibody interactions etc. These interactions 
are expressed as a sum of terms 

^j+ri^j+r2^j+r3 (37) 

where the number of the raising operators should be equal to the num¬ 
ber of the lowering operators from the U(l) symmetry. Then it is easy 
to show the inequality fl35|) . 


3.2 Continuity of energy spectra 


Theorem 3. The lowest energy spectra in the subspace of = {S — 
m/n)L -|- {m and n are coprimes, independent of L; is an 

integer with |AS'|.| L) are continuous as a function of the wavenumber q 
in the inhnite hmit(L —)■ cx)), except = 0. 

Proof. We shall prove for the case ASf = ±1. Generalization is trivial. 


1 . n = 1 case 


We set the left hand side of equation (15^ of theorem 2 as 

5^E{q) = E{Sf] q - 27r/L) + E{Sf] q + 27i/L) - 2E{Sf] q). 

If the energy spectrum were a step function {E{Sf;q) = Q{q 
then from the following relations: 

(0 

1/2 


S E{q) = 


0 


\q - qs\ > 4:Tr/L, 

q = qs- 27r/L, 

q = qs, 


- 1/2 q = qs + 27i/L 


(38) 

qs)), 


(39) 
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the inequality fl5^ could not be satisfied. 

This result denies the step discontinuity. Consequently, possibilities for 
the essential discontinuity and the infinite discontinuity are denied. 

Therefore, the lowest energy spectrum is continuous. 

2. n = 2 case 

Using theorem 2 twice, we obtain 

E{S^- q - Att/L) + E{S^- q + dvr/L) - 2E{S^- q) < 0(1/L), (40) 

therefore, we can prove the continuity similarly to the n = 1 case. 

3. n general 

One can prove the continuity similarly. 


□ 

[Remark] One cannot prove the continuity of the lowest energy spectrum 
in the = {S — m/n)L subspace. 

[Corollary] Although the lowest energy spectra of = ±1, ±2, • • • are 
continuous, the derivative of the spectra may be discontinuous. For example, 
a cusp-like behavior is possible 




9 c) for g > gc, 

Mq-Qc) for g < g^. 


(41) 


in the neighborhood of qc- In this case, from the inequality ([35]), the following 
restriction holds: 


(ui - V2) < 


Q(1) 

27r|S'^ 


(42) 


3.3 Periodicity of energy spectra 

For special values of the magnetization, the wave number change of fl34l) may 
return to the original wave number in finite times. Then using the inequality 
fl32|) several times, we can show the periodicity of the lowest energy spectra. 
This includes the S’!- = 0 case of the original LSM |1] or [m. 
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Theorem 4. (Extension of Oshikawa-Yamanaka-Affleck (OYA) theorem [T3] 

) 

The lowest energy spectra in the snbspace of = {S — m/n)L {m and 
n are coprimes) are periodic with g —)■ g + 27r/n in the inhnite limit: 

|E(A^;g) -E(S'^;g + 27r/n)| < 0{l/L). (43) 

Proof. 1. For the n = 2 case (or S — Sf/L = 1/2, 3/2, • • •) 

From theorem 2, we obtain 

E{Sf; q-7i) + g + vr) - 2E{Sf; q) 

= 2{E{Sf- g + tt) - E{Sf- q)) < 0(1/L), (44) 

where we nsed E(Sf; q' + 2%) = E{Sf; q'). 

Conversely, we can show 

F(5-g)-F(^-g + 7 r)<0 (l/L), (45) 

therefore 

\E{Sf- q) - E{Sf- g + vr)! < 0(1/L). (46) 

2. For the n = 3 case (or S — Sf/L = 1/3, 2/3, • • ■) 

From theorem 2, we obtain 

E{Sf- q + 27r/3) + E{Sf- q - 27r/3) - 2E{Sf] q) < 0(1/L). (47) 

Secondly, applying theorem 2 to the lowest energy state with g + 27r/3, 
we obtain 

E{Sf- q + 47 r/ 3 ) + E{Sf- q) - 2E{Sf- q + 27r/3) 

= E{Sf-, q - 27r/3) + E{Sf; q) - 2E{Sf-, q + 27r/3) < 0(1/L), (48) 

where we nsed E{Sf; q' + 27r) = E{Sf; q'). 

Thirdly, applying theorem 2 to the lowest energy state with g — 27r/3, 
we obtain 

E{Sf] q) + E{Sf] q - 47 r/ 3 ) - 2E{Sf-, q - 27r/3) 

= E{Sf- q) + E{Sf- q + 27r/3) - 2E{Sf- q - 27r/3) < 0(1/L). (49) 
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(50) 


Combining fH7|) x2 + we obtain 

E{S^- q + 27r/3) - E{S^] q) < 0{l/L) 
On the other hand, from (HHj) x2 + (jl9]), we obtain 
E{S^t] Q) - E{S^] q + 27r/3) < 0(1/L) 


3. 


Therefore, 

\E{S^;q) - E{S^;q +27r/3)\ <0(1/1). 
For general case {m,n), we can prove similarly. 


(51) 

(52) 
□ 


Theorem 5. (Second extension of OYA theorem) 

In the inhnite limit, the lowest energy spectra in the subspace = 
L(S — m/n) + AS^ ( AS^ is an integer with |AS'|.| L ) are periodic with 
g —)■ g + 271/n. 

|F;(F^;g) -E(5^;g + 27r/u)| < 0(1/L). (53) 

Proof. 1. For the ASf = 1 and (m,n) = (1,2) case 
From theorem 2 we obtain 

E{Sf-, q-7i + 277/1) + E(Sf\ g + TT - 277/L) - 2E(Sf] q) 

= E(Sf- q + 77 + 277/L) + E{Sf; q + 77 - 277/L) - 2E(Sf; q) < 0(1/L) 

(54) 

where we have used E(Sf;q' + 27r) = E(Sf;q'). In addition, from 
theorem 3, the lowest energy spectrum is a continuous function of g 



\E(Sf; q' ± 277/L) - E(Sf- q') \ < 0(1/L) 

(55) 

Therefore, we 

obtain 



E(Sf- q + 77)- E(Sf- q) < 0(1/L). 

(56) 

Conversely, we can show 



E(Sf- q) - E(Sf- q + 77)< 0(1/L) 

(57) 

In summary 

\E(Sf-q)-E(Sf-q + 77)\<0(l/L). 

(58) 
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(59) 


2. For general AS^ and {m,n) cases, one can prove 


\E{S^- q) - E{S^; q + 27r/n)| < 0(1/L), 


similarly as above and theorem 4. 


[Remark] 



□ 


Figure 1: Elementary excitation of S=l/2 Heisenberg chain 

The periodicity of theorem 4 and 5 is applied only for the lowest energy 
spectrum in not for higher excitations. For example, St = 1 excitation 
for the S=l/2 Heisenberg model, the lower bound behaves as [3] 

J 71 

E = —smq, (60) 

above them there is a continuum of states bounded above by [5] 

E = Jtt sin (61) 

(See Fig. d]) 

3.4 Correlation 

In the case where S is an integer and = 0, one cannot decide whether 
the ground state is gapless or not, with the LSM theorem m- In general 
S — S^/L integer, there is a similar statement [14] . 

However in such a situation, there is a restriction on the expectation value 
of the matrix elements (or correlations). 
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Theorem 6. When = {S — m)L {m: integer), the following ineqnality 
holds. 

(SI; ‘imS)'MS: + - 2H)\SI.; g) > 0. (62) 

Proof. We can show fl62|) . using the inequality fl3^ in theorem 2. 

For example ([T]) with J(r) = this means 

{Sf = {S-m)L;q\S^S^^,\Sf = {S-m)L;q)<0. (63) 

4 Discrete symmetries 

In this section, in addition to the translational symmetry, we will discuss 
discrete symmetries (space inversion, spin reversal). 

4.1 Symmetry in the lowest energy spectrum 

According to the Bethe ansatz for the S=l/2 XXZ spin chain, the lowest 
energy spectrum is symmetric with q = n/2. This can be proven even when 
there is no exact result. 

When the Hamiltonian is invariant under the space inversion (site parity 
or link parity), besides the translational invariance, there is E{Sf]—q) = 
E{Sf] q) symmetry in the energy spectra. 

Theorem 7. (Third extension of OYA theorem) 

In the inhnite limit, the lowest energy spectra of Sf = {S — m/n)L + 
AS"!.; = 0, ±1, ±2, • • •) are symmetric with respect to q = vr/n, 

\E{Sf-, q) - E{Sf- 21,jn - q)\ < 0{1/L). (64) 

Proof. Combining the theorem 5 and E{Sf] —q') = E{Sf]q'), we can prove 
this theorem. □ 

[Remark] 

• From the above theorem, for the S half-integer case, the lowest energy 
spectrum of Sf = 0, ±1, ±2, • • • is symmetric with respect to q = t, 12 
in the inhnite limit 

• This symmetry does not hold for higher energy spectra than the lowest 
spectrum. (See Fig. [1]) 
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4.2 Discrete Symmetries and twisting operator 

4.2.1 Site parity 

Lemma 2. (Space inversion and twisting operators) 

(65) 

Proof. We can obtain 

Pf/S,p = exp ( - ■») 

V i=i 

V j=0 

= U^^^^exp{T27ri{Sf - SL)) exp{T27ri{Sf - S)), (66) 

where we used Sq = Sf. Combining this relation with the fact that the 
eigenvalues of Sf — SL and Sf — S are integers, we can obtain □ 

4.2.2 Spin reversal symmetry 

Lemma 3. (tt rotation around the j/-axis and twisting operator) 

UiU±2. = ( 67 ) 

Proof. 


= exp T— 5^ - S) 

V i=i 

= “p(±¥tw-s)) 

= U%^exp{±2mS{L + l)) 


exp 



( 68 ) 


Combining this relation with the fact SL is an integer, we can show equation 

dnn). □ 
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4.2.3 Link Parity 

Since the twisting operator defined fl2^ and the link parity Punk are not 
compatible with lemma 2, we should define another type twisting unitary 
operator as equation (17) and appendix C in [20] 






Lemma 4. One can show the relation 

p frtwl p _ fjtwl 

-nink(^± 27 r-nink — 


^27r' 


Proof. 


27ri 


= exp Pf-p - l/2)(S£_^+i - S) 


L 

2711 


i=i 

L 


= exp T-p - J + 1 - 1/2)(S| - S) 


i=i 


= P;texp(T27rz(^^-^L)) 


_ -frtwl 

— ^T27r- 


(69) 


(70) 


(71) 

□ 


However, the corresponding relation for lemma 3 should be changed as 

Lemma 5. 

= U%Ul. (72) 

Proof. 

( 2771 ^ 

(UpUtlP = exp U-p 5^(J - l/2)(-S| - S) 


L 

27ii 


1=1 

L 


AttI 


= exp 

'jtwl 


J2u - 1/2)(S' - S) exp ±-f-Sj2(3 - 1/2) 


1=0 

= exp(±27rlS'L) 

_ fjtwl 

~ 12^277 


1 = 1 
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(73) 

□ 







4.3 Combination of discrete symmetries 

Theorem 8. (Corresponding to the eqnation (29) of [TT] i 
On the combination PU'^, the next relation holds 

(PU’WSi, = ( 74 ) 

Proof. Combining lemmas 2 and 3, one can prove this. □ 

Theorem 9. On the combination the next relation holds 

{Pi^^uy)uit = (75) 

Proof. Combining lemmas 4 and 5, one can prove this. □ 


For example, in the spin half-integer case, since the eigenstate l-^o) = 
\Sf = 0;g = 0) satishes PUy\'ipo) = ±|t/’o), therefore from theorem 9 we 
obtain (Pf7^)f/*);'|V’o) = or {^po\Uff\^lJo) = 0. 


5 Consideration of several models 

In this section we consider several models other than ©■ Since the lattice 
strnctnre or the symmetries of them may be different from the model ([T]), we 
should slightly modify the statements of theorems 1-9 in several cases. 


5.1 XXZ spin chain with next-nearest-neighbor inter¬ 
action 

We consider an XXZ spin chain with next-nearest-neighbor (NNN) interac¬ 
tion: 


H = 



+ s-sy.) + AS'sy.) 



+y s;+2)+Ayy+2 

(76) 


This Hamiltonian is invariant with 2 ;-axis rotation, translation, space inver- 
sion(link parity, site parity), and vr rotation around the j/-axis. 

Theorems 1-9 hold. When a > 0 with frustration, the MLM theorem 
does not hold, therefore the uniqueness of the ground state may be broken. 
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In the S=l/2 case, there are several exact results. At a = 0 there is the 
exact results by Bethe ansatz [21 El IH [S] , where energy spectra are known, and 
these results are consistent with our theorems. At a = 1/2 there is another 
type exact result [6], where two-fold exactly degenerate ground states g = 0, vr 
exist even for hnite size. Unfortunately, there is no exact result for excitation 
spectra at a = 1/2. 


5.2 Bilinear-biquadratic (BLBQ) spin chain 

We consider a bilinear-biquadratic (BLBQ) spin chain : 


H = 



(77) 


This Hamiltonian is invariant with z-axis rotation, translation, space inver- 
sion(link parity, site parity), and vr rotation around the y-axis. 

Note that in this case although there appear terms such that 


C+ c— C+ c— C+ c- Q— Q+ QZ QZ Q+ Q— 

W W+IW W+1’ W W+IW W+I’ WW+IW W+I’ 


(78) 


it is straightforward to show the relation fl35p . Therefore, theorems 1-9 hold. 
When a > 0 with frustration, the MLM theorem does not hold. 

In the case of S=l, there are exact solutions by Bethe ansatz at a = 1 
na where gapless excitations are at g = 0,±27r/3, and a = —1 [H] where 
gapless excitations are at g = 0,7r. At a = 1/3 there is another type exact 
result [7|. A unique ground state with an energy gap has been proved. These 
results are consistent with theorems 1-9. 


5.3 XXZ spin chain with staggered field 


" = E + y sy.) + 


i=i 


-lysi 


(79) 


This Hamiltonian is invariant with z-axis rotation, translation by two sites, 
space inversion (site parity). 

Theorems 1-7 hold with several changes. 
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Considering two spins in one nnit cell, we slionld nse the translation 
operator by two sites and rewrite the lemma 1 as: 

i.U,a)^U'Z = y‘:(r>t,i)"exp , (80) 

where we nse 2S as the integer. 

5.4 XXZ spin chain with bond-alternation 

(jl'S/'SJ+i + 'S7'S/+i) + (81) 

This Hamiltonian is invariant with z-axis rotation, translation by two sites, 
space inversion (link parity), and vr rotation aronnd the y-axis. 

Theorems 1-7 and 9 hold with several changes. Lemma 1 becomes 

[UafUif = UZ‘0,„f exp ■ (82) 

In the case of 5 = 0, A 3> 0 the gronnd state is Neel state, that is, two-fold 
degenerate {q = 0 , 7 r). Therefore, when |(5| 1, corresponding to the folding 

of the Brillonin zone, the lowest gronnd state is two-fold degenerate in the 
snbspace = 0; q = 0 bnt different discrete symmetries (Punk = ^ 

versns Punk = — 1 ; = —1 ), which is robnst against pertnrbation 5 7 ^ 0 . 

5.5 Dzyaloshinskii-Moriya type interaction 

The Dzyaloshinskii-Moriya type interaction ns in 

Hdm = X = I E('5/'57+i - Ai) (83) 

j 3 

is U(l) and translational invariant, bnt antisymmetric for the space inversion. 

One can easily show fl3^ . therefore theorems 1-6 hold. Note that in this 
case there is not the q —)■ —q symmetry in the dispersion curve, which may 
be related with the spin spiral ordering. 
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Figure 2: Spin Ladder with nonsymmetric interaction 

5.6 Nonsymmetric spin ladder 

We can consider a spin ladder model with nonsymmetric interactions shown 
in Figl2]with U(l) symmetry. 

One can easily show fl3^ . therefore theorems 1-6 hold. 


6 Spectra minima from the MLM theorem 

For hnite magnetization = {S — m/n)L, {S — m/n)L ± 1, • • • cases, from 
theorems 4,5 and 7 there remain two possibilities for the minimum of en¬ 
ergy spectrum; g = 0 or g = vr. Using only the symmetries such as U(l), 
translation, space inversion or spin reversal, one cannot conclude which is 
appropriate. 

Nevertheless, in the nonfrustrating case, using the MLM theorem [3,19], 
one can distinguish the above two situations. Using the alternating operator: 

(84) 

then we obtain 

(^alt)t^±f;alt _ = S]. (85) 

For example, considering the NNN spin chain fl76|) with a < 0, the unitary 
transformation with this operator is 

L 

(^alt)t^f)alt ^ 

i=i 
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therefore the off-diagonal elements of the Hamiltonian become nonpositive, 
and one can use the Perron-Frobenius theorem (diagonal elements can be 
adjusted by adding a scalar multiple of the identity operator). 

When system size L is hnite, from the Perron-Frobenius theorem, the 
lowest energy state of fl86|) in each subspace is unique with the wavenum¬ 
ber g = 0. Returning to the model fITBD with a < 0 by the unitary operator 
fl84p . the minimum of the energy spectrum is located at g = 0 in the case 
SL — even integer, whereas the minimum of spectrum is at g = vr in the 
case SL — odd integer. Especially when L = 2n, energy spectrum in 
= {S — m/n)L subspace has a minimum at g = 0. Similar discussion can 
be applied for the BLBQ model fl77)) with a < 0. 

Although for frustrating case, the above consideration may become inef¬ 
fective, one may expect such an even-odd difference in the SL — subspace. 

7 Comparison with previous works 

In this section we compare previous works and our theorems. 

7.1 Problem of the original LSM discussion 

For the S=l/2 XXZ model 

" = E (^(4*57+1 + SS+i) + AS'S'+i) , (87) 

Lieb, Schultz and Mattis [1] discussed as follows: 

- H 

L 

= ^ - (^5+4"+i(exp(-^27r/L) - 1) + 4”*5/+i(exp(+z27r/L) - 1)) 
i=i 

= (t) i (»>= (x) -1) 

(88) 
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and 


i=i 


(^0 



,H 


^o) 


0, (89) 


where I'ho) is the ground state. The remaining reasoning is similar to ours. 

However, we think there is a technical problem in their method. For later 
convenience, we rewrite 

-h = ^{^+ a^), 

^ + {UXVHUX - 2 ^, 

AST ^ - iU%jHU%^. (90) 


For flHTl) it is easy to show that the second term of (l90D becomes a com¬ 
mutator form of the Hamiltonian 


AST = i sin 



(91) 


Whereas for the NNN Hamiltonian 0761) . there appears an additive term to 
the commutator: 


AST = i sin 




-|- 2m sin 




(92) 


For the BLBQ model 077)1 . it is cumbersome to derive [TT] 


AST = i sin 


27r 




M 




+ 0(T-^) 


(93) 


These calculations, which become more difficult for complicated models, are 
highly model dependent, which are different from our model independent 
proof of theorem 2 etc. 
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7.2 Problems of Oshikawa-Yamanaka-Affleck discussion 

Oshikawa et al. [H] have assumed the space inversion symmetry or the spin 
reversal symmetry of the model. Although they have not presented explicitly 
the reason, we guess it would be a workaround to avoid the difficulty of the 
previous subsection. In fact, when model is space inversion symmetric, from 
lemma 2 it suffices that 

FASTP = -AST, (94) 

or when model is spin reversal symmetric, from lemma 3 it suffices that 

= -AST, (95) 

therefore, it is apparent that 

(^ol^l^o) = 0, (96) 

where the state Tl'ho) = or = l^o)- 

However, there are problems. First, for the states \S^;q), the discussion 
becomes not so simple; one should use the combination \S^;q) + \S^;—q), 
which leads the three wave numbers inequality fl32|) . not the original LSM 
type. 

Secondly, with the discrete symmetry assumption, several models are 
excluded such as the Dzyaloshinskii-Moriya interaction, or nonsymmetric 
spin ladders, or distorted diamond chain models with a staggered held. 

8 Discussions on energy spectra 

In this section we illustrate possible energy spectra from our theorems. 

8.1 Dispersion curve for frustrating system 

First we treat the frustrate systems with space inversion symmetry. The 
lowest energy spectrum of = ±1 is continuous from theorem 3, and there 
are symmetry restrictions from theorems 1,5 and 7. 

• S: half-integer spin 

Regarding the dispersion curve in = ±1, there is the possibility of 
the four lowest points, other than the conventional two lowest points 
q = 0,71. (See hgure E]) 
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Figure 3: Possible energy spectrum S = 1/2, 3/2,--- (a) two lowest points 
(b) four lowest points 


For example, when the next-nearest-neighbor interaction a is large 
enough in the S=l/2 NNN XXZ spin chain fl76l) . dispersion curve may 
have four minima. 

• S: integer spin 

Regarding the dispersion curve in = ±1, there is the possibility of 
the two lowest points(g 7 ^ 0 , tt) , other than the conventional unique 
lowest point. (See figure 0]) 


-n q=0 n -n q=0 n 

Figure 4: Possible energy spectrum for S' = 1, 2, • • • (a) unique lowest point 
(b) two lowest points 

For example, when the biquadratic interaction a is large enough in the 
S=1 BLBQ spin chain (1771) . dispersion curve may have two minima. 

8.2 Magnetization plateaux 

We consider the energy spectra at magnetization plateaux. As examples, we 
treat the = {S — 1/3)L case and the = {S — 1/3)L -|- 1 case. In these 
cases, the lowest energy spectra are periodic with the wavenumber q ^ q + 
27r/3. In the = {S — 1/3)L subspace, the lowest energy spectrum may be 
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discontinuous, and the lowest energy state is located at g = 0 from the MLM 
discussion in the section 6 (£gure[S](a)). In the = (S' —1/3)L + 1 subspace, 
the lowest energy spectrum must be continuous with q from theorem 3, and 
the minimum point is at g = vr (£gure[5] (b)). 



Figure 5: Possible energy spectra (a) = {S—1/3)L (b) = (S —1/3)L+1 


8.3 Nonsymmetric dispersion 

We consider energy spectra of a model with Dzyaloshinskii-Moriya type in¬ 
teractions flH^ or the nonsymmetric spin ladder (see figure [2]). In this case, 
energy spectra are periodic from theorem 4 and 5, whereas they are not 
symmetric under g -H- — g (hgureE]). 



Figure 6: Possible energy spectrum for nonsymmetric model 


9 Conclusions 

We have extended the LSM theorem including the frustrated case, because 
we have not used the uniqueness condition of the lowest state in each 
subspace. We have also extended the LSM theorem for nonsymmetric case, 
for example, the Dzyaloshinskii-Moriya interaction. 
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Regarding the continuity of the lowest energy spectra in the = ±1, ±2, • • • 
subspace, we have completed insufficient points of previous works [13], since 
they used inequality relation only once. Relating to this, Oshikawa [21] stated 
“that the incommensurate hlling gives a gapless spectrum is empirically rec¬ 
ognized more or less”, and he tried to prove it with topological arguments. 
Although his argument is intuitive, it is not convincing. 

There are several remarks. First, on the ordering of the lowest energies 
of each subspace; one can say nothing about it from the extended LSM 
theorem itself. Secondly, the lowest energy states in the = {S — m/n)L 
subspace may be discontinuous as a function of wave number q, which is 
natural, considering Neel or dimer states or magnetic plateaux. Thirdly, 
cautions should be taken on the number of spins in the unit cell. 

Finally, the original LSM theorem has also been extended for fermion 
systems on the lattice [22l |23] . It will be interesting to consider our methods 
for fermion models with frustration. 
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